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Abstract 

We present a simple method for calculating the energies and the widths of quasiparticle resonant 
states. The method is based on BCS equations solved in the Berggren representation. In this 
representation the quasiparticle resonances are associated to the Gamow states of the mean field. 
The method is illustrated for the case of neutron-rich nuclei 20 ~ 22 O and 84 Ni. It is shown that the 
contribution of the continuum coupling to the pairing correlations is small and largely dominated 
by a few resonant states close to the continuum threshold. 
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I. INTRODUCTION 



In nuclei close to the drip lines the low-lying excited states belong mostly to the continuum 
part of the spectrum. It is therefore not surprising that a lot of effort has been devoted 
recently to the treatment of the continuum in nuclear models. 

One of the first attemps to evaluate continuum configurations in nuclei was done within 
the continuum shell-model ( see Ref.jjJ and the references therein). In this shell model 
approach the coupling to the single-particle continuum is introduced through the scattering 
states of real energies. A different approach for evaluating the continuum processes in the 
shell-model framework is based on Berggren representation. In this representation the single- 
particle continuum is formed by a finite number of Gamow resonances and a continous set 
of scattering states of complex energy. The advantage of a shell model based on Berggren 
representation is that it allows a direct calculation of multi-particle resonant states in nuclear 
systems Q, S, fl E , [J . 

The continuum part of the excitation spectrum was also investigated using mean field 
techniques. One of the most known of these approaches is the continuum- RPA [?], Is| . In 
this approach, well suited for the description of particle-hole (p — h) excitations in closed 
shell nuclei, the contribution of the continuum to the p — h response function is calculated 
exactly. To single out from the continuum p — h response the states corresponding to the 
p — h resonances, a resonant-RPA approach was also developed 0, ^| , in which the coupling 
to the single-particle continuum is introduced through a finite set of Gamow functions. 

For even-even nuclei with open shells, where pairing correlations are important, the con- 
tinuum part of the spectrum corresponds to two-quasiparticle excitations. The unbound 
two-quasiparticle excitations can be described in the continuum-QRPA approach 
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Self-consistent continuum-QRPA calculations are presently done only with zero-range forces 



12j . At variance with the continuum- RPA, in these calculations the continuum contribution 
is taken into account only up to an energy cut-off. 

In odd-even nuclei the most simple excitations are of the one quasiparticle type. The 
continuum part of the one quasiparticle spectrum can be calculated by solving the HFB 
equations with scattering type boundary conditions |l3[- By using the behaviour of the 
phase shift close to a resonance, one can find the quasiparticle resonances of physical interest. 
The quasiparticle resonances can be found also by solving the HFB equations in complex 
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energy plane, as have been already done for some particular hamiltonians jl^j ]. 

However, the simplest framework to evaluate one quasiparticle resonances is the BCS 
approach. The contribution of resonant states to pairing correlations can be introduced 
in the BCS equations through the scattering states with energies located in the vicinity of 
single-particle resonance energies 16, ^7|. Alternatively, the resonant continuum was 
studied by solving the BCS equations in a basis containing only bound and Gamow states 



18l 19^ . The drawback of this study is that the global quantities associated to the bound 
state of a nucleus ( e.g., binding energy, radii) become complex. The scope of this paper 
is to show how one can calculate within the BCS approach the energies and the widths 
of quasiparticle resonances by taking into account the whole contribution of single-particle 
continuum to pairing correlations and, in the same time, by preserving the correct properties 
of physical observables. 

II. FORMALISM 

A. Berggren representation 

The quasiparticle spectrum is described here within the Berggren representation. Since 
this representation is not very common in nuclear structure calculations, we will first describe 
it briefly, emphasizing those features that we will need in the present paper. For more details 
see Refs. flHQ- 

The Berggren representation is formed by bound states, Gamow resonances, and a con- 
tinuous set of scattering states of complex energy. For a given mean field potential (e.g., of 
Woods-Saxon type), the Gamow states are the outgoing solutions of the Schrodinger equa- 
tion and correspond to a discrete set of complex wave numbers k u = k v — i^ u . To each 
Gamow state u v ij(r) with the wave number k v — k v — ij u is associated a "mirror" state 
Uuijij) with the wave number k u = —k v — i^ V) such that w*^-(r) = u u ij{r). The Gamow 
states and their mirrors form a biorthogonal set which can be normalized to unity by using 
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various regularization schemes 



23j . All the wave functions written above are the radial 
single-particle wave functions of angular momentum Ij. 

The scattering states introduced in the Berggren representation belongs to a contour in 
the complex momentum or complex energy plane. An example of such a contour in the 
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FIG. 1: Contour on the complex energy plane representing the continuum path. The points V{ are 
the vertices defining the contour. The open circles labelled by Gi indicate the complex energy of 
the Gamow resonances enclosed by the contour while the open squares label by Bi are the bound 
states. The dot is a pole of the S-matrix excluded from the Berggren expansion 



complex energy plane is shown in Figure ^ 

The contour encloses the poles associated to the Gamow states which are taken explicitely 
into the Berggren representation. The scattering states belonging to the contour, uij(e,r), 
the enclosed Gamow resonances, u u ij(r), and the bound states, u n ij(r), form a complete set, 
i.e., 

5(r - r ') = u n ij(r) u n ij(r') + ^ u„ij(r) w^(r') + de uij(e,r) utj(e,r'). (1) 

n v L 

In numerical applications the integral along the path L is evaluated by using a finite number 
of scattering states of energy E{ defined by a chosen quadrature rule. In what follows this 
set of scattering wave functions together with the bound and Gamow states are denoted by 
a unique function, i.e., 

{Uuj(r) bound states and resonances 

y/XiUij(Si,r) scattering states , 

where the quantities Xi depends on the quadrature. In general, they are defined by Xi — ^iLi, 
where uii are the quadrature weights and Lj are the derivative of the complex contour with 
respect to the parametrisation variable. 

In principle, a Berggren representation can be formed by using any set of Gamow states 
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enclosed by a contour in the complex plane. However, from physical point of view it is better 
to choose a basis which includes explicitely only narrow Gamow resonances. The contribu- 
tion of wide resonances which are not included explicitetly in the Berggren representation 
is automatically taken into account through the scattering states belonging to the contour. 



B. BCS equations in Berggren representation 

The BCS equations in Berggren basis will be derived here from the Hartree-Fock- 
Bogolibov (HFB) equations in coordinate representation. For zero range forces and for 
spherically symmetric systems the radial HFB equations have the following form j^J: 

h(r) - A A(r) \ / U k (r) \ = I U h (r) \ 
A(r) -h(r) + A J \ V k (r) ) * \ V k (r) ) ' 

where A is the chemical potential, h(r) and A(r) are the mean field hamiltonian and pairing 
field, respectively. U k {r) and V k (r) are the up and down components of the HFB wave 
function, denoted below by >). At large distances U k (r) decays exponentially for E k < 
— A and behaves as a scattering state for E k > —A. On the other hand, the function 
Vfc(r) goes to zero exponentially for any value of the quasiparticle e nerg y (for a general 
discussion on the asymtotic properties of the HFB wave functions see |25|). Consequently, 
the quasiparticle spectrum is discrete for E k < —A and continous for Ej. > —A. In the 
ground state of the system the particle density and the pairing density are given by: 

-°( r ) = i^E( 2 ^ + 1 )^( r )^( r ) • ( 3 ) 

4 7T z — ' 

k 



p{r) = ^Y, {2 ^ + l ) u ^)y^) ■ ( 4 ) 

4 7T z — ' 

k 

The summations in this equation is over the quasiparticle spectrum of the system (for 
E k > — A the sum should be replaced by an integral) up to a given cut-off energy. 

In the calculations presented in the next section we take for the mean field a Wood-Saxon 
potential. The pai ring field is calculated self-consistently by using a pairing interaction of 
the following form 26]: 

V(t - r') = V [l - v(-) a \S(r - r') = V eff 5(r - r'). (5) 

Po 
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With this pairing interaction the pairing field is given by: 

A(r) = ¥fp{r). (6) 

Before going further with the derivation, it is worth mentioning how the particle and 
pairing densities behave in excited states belonging to the continuum part of the quasiparticle 
sectrum of odd-even nuclei. Formally, in order to get the excited state coresponding to a state 
k of an odd-even nucleus one can simply interchange the eigenvector (Uk,Vk) by (V£,U}1) 
Consequently, the particle density corresponding to an excited state k will depend 
on \Uk\ 2 . Thus, if the excitated state belongs to the continuum, the function XJ k is of a 
scattering type; hence, there is a finite probability for the particle to be found at infinity, 
as physically required. On the other hand, by interchanging U). and one does not affect 
the asymtotic behaviour of pairing density, which is still going to zero exponentially at very 
large distances. This is again a physically meaningful condition since when a particle is very 
far from the nucleus, its contribution to the pairing correlations should vanish. 

Starting from the HFB equations in coordinate space we can easily get the HFB equations 
in Berggren representation. In order to accomplish that, we first expand the HFB wave 
functions in the Berggren basis: 

U k (r) = J2< ^\ Uk > ^( r ) = 2 U ^ k ^ ( ? ) 

i i 
i i 

where the expansion coeficients, related to the Bogoliubov transformation, are complex 
quantities. Then one multiplies the HFB equations from the left with the complex conjugate 
of a mirror Berggren vector, ip*(r) and one performs the radial integration. One thus gets 
the HFB equations in Berggren representation: 

, A A fJ - \ ( Uu \ / U * 
X.j -Si + XJ \ V id ) 1 \ Vij 
which formally have the same structure as in any hermitic representation. The essential dif- 
ference is that the matrix elements of the pairing field are defined with the Berggren metric, 
i.e., by using mirror states in the bra positions and employing regularisation techniques for 
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calculating the diverging integrals. Thus, Ay is defined by: 

A itj =< ^|A(r)|V>j >= J drr 2 ^(r)A(r)^-(r) , (10) 

where tpi is the mirror state corresponding to fy. 

Finally, from Eq. Q one can get the BCS equations in the Berggren representation 
by neglecting the off-diagonal matrix elements of the pairing gap matrix. Physically, this 
approximation means that one neglects the pairing correlations associated to the Cooper 
pairs formed in states which are not time-reversed. In this approximation the expansions 
(7,8) of the HFB wave function are reduced to one term. The corresponding expansion 
coeficients, denoted below by and Vi, are complex quantities. They satisfy the condition 
u? +v? — 1, which is obtained from the Berggren normalization condition of the HFB wave 
function i.e., < ^il^j >= 1. By combining this condition with Eq.fjOJ one gets the BCS 
equations in the Berggren representation: 



1 q £i — A, 



E; 



>< ~ A) 2 + A 2 



where 



A; = J dr r 2 ^ 2 (r)A(r) 
and the chemical potential is found from the particle number equation: 



N 



(2* + 1) v} tf(r) 



r dr 



(11) 
(12) 

(13) 
(14) 



In conclusion, apart from the radial integrals which are calculated with the Berggren 
metric (and a given regularisation procedure), the BCS equations in the Berggren represen- 
tation have formally the same expression as in any hermitic and discrete basis. However, 
one should keep in mind that the pairing gaps and the occupation probabilities associated 
to the scattering states correspond to an energy interval. 

The advantage of working in Berggren representation is that the quasiparticle resonances 
appear as unique states, associated to the Gamow functions. In this way the widths of 
the quasiparticle resonances, given by the imaginary parts of the quasiparticle energies, are 
calculated unambiguously. 



In principle, one can get the complex energies associated with the resonant states by 
solving the BCS equations in complex energy plane. However, due to the strong non- 
linearity of the BCS equations, this is a difficult numerical task. One can avoid this task 
by using the fact that the pairing field and the chemical potential do not depend on the 
representation. Thus, they can be calculated by solving the BCS equations in a real energy 
representation, which can be obtained by choosing the L contour to be the real energy 
axis. From the pairing field calculated using the real energy axis, we can then get, by using 
Eq.(llc), the pairing gaps associated to the Gamow states. With these values of the pairing 
gaps, which are complex quantities, one can finally calculate, by using Eq.(llb), the energies 
of the resonant quasiparticle states. 

III. APPLICATIONS 

The scope of the calculations presented below is twofold: a) to illustrate, within the 
approximation discussed above, how the energies and the widths of quasiparticle resonant 
states behaves in neutron-rich nuclei; b) to analyse what is the relative contribution of res- 
onant and non-resonant continuum upon pairing correlations in nuclei close to the dripline. 
For these purposes we have performed calculations for the isotopes 20_22 O and 84 Ni. 

A. Mean fields and single-particle states 

The BCS calculations presented in this section are based on a mean field of Woods- 
Saxon form. In addition to the mean field, we take also a standard spin-orbit interaction, 
with the form factor given by the derivative of a Wood-Saxon function. The Wood-Saxon 
parameters we used for oxygen isotopes have the following values: Vo=55.8 MeV, V^ Q =12.12 
MeV, r =1.21 fm and a = 0.65 fm. The values of r and a are the same for the mean field 
and for the spin-orbit interaction. With these parameters we get for the isotope 17 the 
single-particle energies shown in Table HJ 

One can notice that the energy of the single-particle state is very close to the exper- 
imental value, equal to (0.94 — i0.05) MeV |28( . The relevant single-particle resonances in 
oxygen isotopes are the states 2d^/2 and 1/7/2- The corresponding complex energies shown 
in Table|Uare the energies of the Gamow states, i.e., e u = e u — iT u /2, where T u is the width 
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TABLE I: Single- particle energies (in MeV) corresponding to 17 and 79 Ni 



State 


17 


— J±.±\->±. £^J.\^kJ 




79 Ni 


Energies 


I/7/2 




7.45-i 1.53 


l/ill/2 




3.56-i 0.02 


1^3/2 




0.91-i 0.05 


157/2 




1.93-i 0.01 


2si/ 2 




-3.26 


2c?3/2 




0.51-i 0.06 


14/2 




-4.12 


3Si/ 2 




-0.82 


lPl/2 




-14.57 


24/2 




-1.44 



of the resonance. 

For the isotope 84 Ni we have chosen the following Wood-Saxon parameters: Vo=42.0 
MeV, V SO =1Q.O MeV, r =1.27 fm and a = 0.67 fm. As in the previous case, the values of Tq 
and a are the same for the mean field and for the spin-orbit interaction. The single-particle 
states corresponding to the major shell N=50-82 are given in Table It can be seen that at 
variance to the situation in stable nuclei ( e.g., tin isotopes) the states ld 3 /2, 1*77/2 and lh u / 2 
are not bound states but single-particle resonances. In the presence of pairing correlations 
these states will generate the quasiparticle resonances discussed in the subsection C. 

B. The BCS solution: contribution of resonant continuum 

With the single-particle states described in the previous section we have solved the BCS 
equations for the pairing force given by Eq.(5). For the parameters of the pairing force 
we have used the following values: a) oxygen isotopes: Vo=-456 MeV fm -3 , 77=1, a=l; b) 
nickel isotope: Vo=-1130 MeV fm -3 , 77=1, a=l.In both cases we have introduced in the 
calculations all the bound states. In order to analyse various approximations of continuum 
treatment, we have performed two types of BCS calculations. 

First, we have considered in the BCS equations the whole contribution of the continuum 
states up to a cut-off energy equal to 10 MeV. The continuum contribution is introduced 
through the scattering states of real energy. The number of scattering states per energy 
interval was increased, especially in the region of single-particle resonances, up to the con- 
vergence. In the calculations we have considered for the continuum states the same (lj) 
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FIG. 2: The scattering state ^3/2 evaluated at the position of the ld^n resonance, i. e. 0.910 MeV 
(full line). The dashed line shows the real part of the corresponding Gamow function, with the 
energy equal to (0.910,-0.050) MeV. 



TABLE II: Average gap, pairing energies (in MeV) and root-mean-square radii (in fm) calculated 
by using the representations discussed in the tex. 







cBCS 






rBCS 




bBCS 








< A > 


E p 


< r > 


< A > 


E p 


< r > 


< A > 


E p 


< r > 


20Q 


1.831 


-1.044 


2.916 


1.800 


-1.013 


2.913 


1.784 


-1.002 


2.912 


22 Q 


1.412 


-0.876 


3.040 


1.384 


-0.856 


3.038 


1.370 


-0.849 


3.038 


84 M 


1.493 


-0.644 


4.650 


1.390 


-0.610 


4.647 


1.386 


-0.606 


4.646 



values as for the bound states. These calculations are called below continuum-BCS (cBCS). 

In the second BCS calculation we have kept from the continuum only the contribution 
corresponding to the resonant states. More precisely, we have included in the calculations 
only the scattering states with the energies in the interval e u ±2T U , where e v and Y v are the 
energies and the widths of the resonances given in Table |l] These calculations are refered 
to as resonant-BCS (rBCS) lrjj ]. 

The results of these two calculations are shown in Table |Hj In order to see what is the 
total contribution of the continuum, in Table |H] we give also the results obtained when in 
the calculations we include only the bound states (bBCS). 

From Table ILT1 one can first notice that for all the isotopes analysed here the contribution 
of the continuum to the correlation energies is rather small. Secondly, one can see that 
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the contribution of the continuum is essentially given by a few low-lying single-particle 
resonances. In the presence of pairing correlations they generate the resonant quasiparticles 
states, analysed in the next subsection. 

C. Quasiparticle resonances 

The quasiparticle resonances, which describe unbound excited states in odd-even nuclei, 
are obtained from the BCS equations written in Berggren representation. In this repre- 
sentation each quasiparticle resonance is described by a unique state, corresponding to a 
single-particle Gamow state. This is different from the case of real energy representations, 
where a quasiparticle resonance is charaterized by a continuous set of scattering states with 
energies in the vicinity of the resonance energy. 

As it was already discussed in the previous section, in order to get the energies and the 
widths of the quasiparticle resonances is necessary to calculate the pairing field and the 
chemical potential. Since these two quantities do not depends on the representation, we can 
calculate them by solving the BCS equations by using a real energy representation. Thus, 
we have solved the BCS equations by using all the bound states and the whole contribution 
of the continuum states up to a cut-off energy equal to 10 MeV. The continuum contribution 
was introduced through the scattering states of real energy. This calculation corresponds to 
the case cBCS described in the previous subsection. 

With the mean field and the chemical potentials provided by the cBCS equations we 
then calculated, by using Eqs.(llc), the pairing gaps associated to a Gamow state and the 
corresponding quasiparticle energies and widths. The results are shown in Table IIII1 From 
these results one can see clearly how the widths of single- particle resonances are changing 
by the pairing correlations. This direct accesss to the widths of quasiparticle resonances 
is the main advantage of expressing the BCS solution in terms of single-particle Gamow 
functions. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we have discussed how the energies and the widths of quasiparticle reso- 
nances can be calculated by using the BCS equations written in Berggren representation. 
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TABLE III: The quasiparticle energies, Ei, the corresponding single-particle energies measured from the chemical potential, i.e., ei = \ei — A| 



and the occupation probabilities. The energies are in MeV. 



State i 




20 






22 




State i 




84 Ni 








Et 


»? 


e; 


Ei 


«? 




e« 


Ei 


«? 


1^5/2 


(0.301,0) 


(1.656,0) 


(0.591,0) 


(1.060,0) 


(1.657,0) 


(0.820,0) 


24/2 


(1.255,0) 


(1.720,0) 


(0.864,0) 


2si/ 2 


(0.564,0) 


(1.301,0) 


(0.283,0) 


(0.195,0) 


(1.020,0) 


(0.596,0) 


3si/ 2 


(0.631,0) 


(1.044,0) 


(0.801,0) 


1^3/2 


(4.736,-0.141) 


(4.878,-0.084) 


(0.015,-0.003) 


(3.977,-0.148) 


(4.082,-0.075) 


(0.013,-0.003) 


24/2 


(0.701,-0.050) 


(1.050,-0.183) 


(0.171,-0.030) 


!/7/2 


(11.273,-1.528) 


(11.310,-1.574) 


(0.002,-0.001) 


(10.514,-1.528) 


(10.541,-1.559) 


(0.001,-0.001) 


!/7/2 


(2.121,-0.010) 


(2.486,-0.021) 


(0.073,-0.002) 
















1^11/2 


(3.745,-0.018) 


(3.937,-0.025) 


(0.024,-0.001) 



In this representation the pairing gaps associated to the quasiparticle resonances, A„, are 
obtained by integrating the pairing field with the square of Gamow states. On the other 
hand, the pairing field, which does not depend on representation, is calculated by solving the 
BCS equations in a single-particle basis formed from bound states and scattering states of 
real energies. Then, the energies associated to the quasiparticle resonances E v are calculated 
by using the standard BCS formula,i.e., E v = y/ {e v — A) 2 + A 2 , where e v are the energies of 
the Gamow states. Since both e„ and Aj, are complex quantities, the quasiparticle energies 
are also complex. As in the case of single-particle resonances, the imaginary parts of E v are 
associated to the widths of the quasiparticle resonances. 

The calculation of quasiparticle resonances was illustrated for the case of the isotopes 
20-22Q an( j 84j^ W as shown that in these isotopes the contribution of the single-particle 
resonant states upon pairing correlations is dominant compared to the contribution of back- 
ground continuum. The most relevant aspect of the BCS calculations presented here is that 
they describe the excitations in resonant quasiparticle states in term of individual states, as 
in the case of bound states excitations. 
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